Quantization of exciton in magnetic field background 
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The possible mismatch between the theoretical and experimental absorption of the edge peaks 
in semiconductors in a magnetic field background may arise due to the approximation scheme used 
to analytically calculate the absorption coefficient. As a possible remedy we suggest to consider 
nontrivial boundary conditions on x-y plane by in-equivalently quantizing the exciton in background 
magnetic field. This inequivalent quantization is based on von Neumann's method of self-adjoint 
extension, which is characterized by a parameter E. We obtain bound state solution and scattering 
state solution, which in general depend upon the self-adjoint extension parameter E. The parameter 
E can be used to fine tune the optical absorption coefficient -R"(E) to match with the experiment. 

PACS numbers: 03.65.Ge, 71.35.-y, 71.35.Cc, 03.65.Db 



I. INTRODUCTION 

Exciton is a bound pair between an electron and a 
hole in a crystal, which interact through Coulomb attrac- 
tion between them. When an electron from the valence 
band is excited into the conduction band, the missing 
electron in the valence band leaves a hole with opposite 
electric charge. The exciton results from the binding of 
the electron with its hole. As a result the exciton has 
slightly less energy than the unbound electron and hole. 
It has an important role to play in the fine structure of the 
absorption edge in semiconductor. It has been seen that 
as the photon energy goes above the band gap energy 
of the semiconductor, the absorption coefficient shows 
spikes due to the formation of exciton [2]. This is an 
important study no doubt, because it explores the band 
structure and the character of the one electron states of 
the semiconductor. The properties of a semiconductor 
are mostly determined by the electronic states which lie 
just below and just above the forbidden gap. 

An appreciable amount of studies have been made to 
explore the band structure of semiconductors by examin- 
ing direct and indirect optical absorption in semiconduc- 
tors Q . This absorption may arise due to the transition 
between the states of the same side of the forbidden gap 
or across the forbidden gap. 

There have been a lot of studies in this field both theo- 
retically and experimentally, see for example [H, [5| . Since 
the constituents of exciton are electrically charged, the 
magnetic fields have certain effects on the energy spec- 
tra of the exciton and it in turn affects the fine structure 
of the absorption. The magnetic fields have a Zeeman 
like effect on the energy spectra of an exciton, which 
has been investigated by studying optical absorption in 
many substances, for example, in InSb Ge and InAs. 
Magneto-optic effect in strained and unstrained Ge have 
been performed in [t']. In [s^, a theoretical study of the 
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effect of magnetic field on the absorption edge in solid 
has been done, neglecting the Coulomb interaction be- 
tween the electron and the hole. It is shown that a se- 
ries of peaks may form due to the transition between 
the sub-bands, which originate due to the splitting of 
electron energy states in the solid by the magnetic field. 
Instead of forming peaks, in some cases a series of steps 
are formed due to magnetic field. But Coulomb interac- 
tion certainly has effect [9j on the absorption in solids. 
It is therefore required to consider Coulomb interaction 
and magnetic field simultaneously in the study of ab- 
sorption In general this problem is hard to tackle 
analytically for any magnetic field. But for large mag- 
netic field an approximation scheme (lo| can be used to 
handle the situation and can be checked experimentally 
for some substances like InSb, Ge etc. According to this 
approximation scheme, which is valid in large magnetic 
field, the Coulomb term should only affect the dynamics 
of the exciton in the direction of the applied magnetic 
field and it is therefore possible to reduce the time inde- 
pendent three dimensional Schrodinger equation into one 
dimensional eigenvalue equation in the direction of the 
magnetic field. The inclusion of the Coulomb interaction 
shows that the the absorption peaks in each magnetic 
sub-band transition is moved towards the lowest exciton 
line. The other effects are the reduction of intensity of the 
continuous absorption and to smooth out the peaks near 
the absorption edge. The peaks are the exciton peaks 
instead of the magneto-optic peaks. 

It is expected that the analytical theory of absorption 
in solids will not be able to exactly explore the true 
nature of the band structure due to the approximation 
used in the calculation. In this article, we therefore try to 
incorporate some short range interactions into the bound- 
ary condition, imposed on the domain of the system so 
that the results based on the approximate model agrees 
with the experiments as much as possible. This nontriv- 
ial boundary condition will in general change the energy 
spectrum of the system and it will depend on a param- 
eter, called self-adjoint extension parameter. The fine 
structure of the optical absorption edge in a crystal due 
to exciton essentially depends on the energy spectrum 
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of the exciton. It is therefore expected to change if the 
spectrum changes. To get a consistent boundary con- 
dition and the domain, so that the Hamiltonian remains 
self-adjoint, we perform a self-adjoint extension (SAE) of 
the one dimensional Hamiltonian in the direction of the 
magnetic field by using von Neumann's method . Self- 
adjointness of a Hamiltonian is essential, because other- 
wise the Hamiltonian would generate complex eigenval- 
ues and the time evolution of the states will not be uni- 
tary. SAE has received lot of interests in recent years 
and is now being used extensively in different branches 
of physics, to explore the nontrivial quantum behavior of 
different systems [H [H [11, [H, [11, E3, [l|] . 

The present article has been organized as follows: In 
sec. II, we discuss about the exciton in a magnetic field 
background and write an effective radial eigenvalue equa- 
tion following '10'], which we need to quantize using non- 
trivial boundary condition. The method of self-adjoint 
extensions (SAE) is discussed in Sec. HI, which is essen- 
tial to get nontrivial boundary condition for the effective 
Hamiltonian H^^^ so that the Hamiltonian is self-adjoint. 
The bound state and scattering state solutions are dis- 
cussed in Sec. IV and V respectively. We conclude the 
paper with a discussion in Sec VI. 



II. EXCITON IN MAGNETIC FIELD 
BACKGROUND 

We consider dynamics of electron and hole pair with 
effective masses rrie and rrih respectively in a medium of 
dielectric constant k. In a magnetic field background B 
the time independent Schrodinger equation for the sys- 
tem is given by [lOj 
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where A{re) and A{rh) are magnetic vector potentials 
at electron and hole positions respectively. 

Considering the fact that the states of interest in opti- 
cal transitions are only those in which k = ke + kh = 
(where fee and kh are electron and hole wave vectors 
respectively). Lamb [l9j simplified the above equation 
by changing the variables to the center of mass coor- 
dinate p — mf.rf. + nihrii/m^ -f m/j and relative posi- 
tion r — re — Vh and using the substitution ^'(re,r/i) — 
[/(r)exp[~^B X p.r] 
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where n is the reduced mass n = iriemh/ {me + nih) of the 
electron-hole system. The presence of Coulomb interac- 



tion spoils the separability of ([2]) in any co-ordinates. So 
it is difficult to solve the equation in any coordinate sys- 
tem. It can however be reduced in an one dimensional 
eigenvalue equation in large magnetic field limit. Note 
that the magnetic field affects the motion of a charged 
particle in a direction perpendicular to the magnetic field. 
In our case we take the magnetic field in z direction. So, 
in z direction the charged particle motion will be a free 
particle motion in absence of Coulomb force. But in pres- 
ence of Coulomb force the free motion in z direction will 
be modified according to Ref. The potential along 

z direction will now be the average Coulomb potential 
over the x-y plane at each point in z axis 10]. From now 



on, we take ti^ — 2fi 



1 , for the simplicity of calcu- 



lations. In these units, the one dimensional Differential 
Equation is of the form [l3| 
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where a and A are constant parameters (see Ref. [ifll for 
details) and its values will depend on the eigenfunctions 
on the perpendicular x-y plane. We thus lead to a one di- 
mensional differential eigenvalue problem which depends 
on the modulus of Cartesian co-ordinate z due to the re- 
fiection symmetry on the z axis. The inner product is 
now defined as 

/"OO 

(Xi(|z|),X2(k|))= / Xi{\z\)x2i\z\)d\z\ (4) 
Jo 

where d\z\ is the measure. In order to find the bound 
state eigenvalue and corresponding eigenfunction of ([3]), 
we need to consider a physically meaningful boundary 
condition. We want to find out this boundary condi- 
tion by von Neumann's method of self-adjoint extension 
(SAE). 



III. VON NEUMANN'S METHOD OF SAE 

Before going into the actual Hamiltonian H^z\ of ([3]), 
it is essential to discuss about the general properties and 
method of self-adjoint extensions of an operator B defined 
over the elements of a Hilbert space TC. Let us consider 
'D{B) as the domain of B such that the operator B be- 
comes symmetric in that domain. The operator B is said 
to be symmetric (or hermitian) if it satisfies the condition 
(^01, 02) = (0i,S02), V0i,02 € ViB), where (., .) is the 
inner product defined over the Hilbert space. Given the 
operator B and the inner product (., .), one can calculate 
the adjoint operator B* from the Green's formula [2l] |. 
This adjoint operator B* should have a domain P(B*), 
which can be obtained from (B* (f>i , (f>2) = i4>ii I3(f)2), 
\f(f>2 e V{B), where (^i € V{B*). Now, if the operator is 
self-adjoint then V{B) = V{B*) and for non self-adjoint 
case V{B) ^ ViB*). For non self-adjoint operators we 
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FIG. 1: (color online) A plot of the function of the eigenvalue 
as a function of the parameter a. This type of plot is 
previously obtained in Ref. 13, 23]. This discussion and the 
properties of the function can be found in [23 |. 



where a± = ±^ie^'''^^, m = {\ + 4r) and W is the 

Whittaker's function ^2?|. The square integrability of 
([5|) can be easily checked from the asymptotic and short 
distance behavior. Note however that the potential at 
\z\ — > is always finite so there is no problem of short 
distance square-integrability but the short distance be- 
havior will help to get the spectrum for the problem. 
The asymptotic behavior is given by 



lim (t>±{\z\) 

|2|— i-OO 
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-(a+|z|)),(6) 



which goes to zero asymptotically, i.e., (/'id^D — > 0. In 
order to get short distance behavior we need to express 
the Whittaker's function in terms of Kummer function 
as 



need to find out self-adjoint extensions, because accord- 
ing to Stone's theorem self-adjointness of an operator 
is necessary in order to guarantee unitary evolution. We 
follow von Neumann's method for finding self-adjoint ex- 
tensions of the Hamiltonian, if there is any. According 
to this method, we first need to find out the solutions 
of the equation B*(jr^ = ±i(jr^. (fr^ are called deficiency 
space solutions, because the numbers are a measure 
by which a non self-adjoint operator is away from self- 
adjointness, where iv^ denotes the number of (fr^ solu- 
tions. For n"*" = n~ = 0, the operator B* is self-adjoint. 
For ri+ = = n ^ 0, the operator B* is not self-adjoint, 
but admits self-adjoint extensions characterized by "n? pa- 
rameters of U (n) . Finally, for ?t,+ ^ n~ , the operator B* 
is not self-adjoint and can't be made self-adjoint. The 
self-adjoint extension of the operator B* is given by the 
domain Vu(n){B) ee <l) + <j)+ + U{n)<j)- , where </> £ V{B). 
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where M is the Kummer function [24]. The short dis- 
tance behavior can now be evaluated as 



lim 0±(U| 
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where A± — 



(l+2m)r(l/2 — a-)- — m)r( l+27n) 



and B± 



n TT (l-|-2m)r(l/2-a±+m)r(l-2m) ' 

Since in our case n+ ~ n_ = 



1, is not self-adjoint 
but it admits self-adjoint extensions, characterized by a 
parameter E and the domain over which the Hamiltonian 
is self- adjoint is given by 



IV. BOUND STATE SOLUTION OF EXCITON 

We now come to the discussion of the Hamiltonian ([3]) , 
for which we seek a self-adjoint domain. We construct an 
initial domain (p{\z\) G where (pil^l) a-^d 0'(|2:|) 

are absolutely continuous and (/)(0) — 0'(O) = 0. It is 
implied that (f>{\z\) and i?|2|0(|z|) are square-integrable 
functions, i.e., belong to the Hilbert space. Note that 
2?(if|2|) is so restricted that the Hamiltonian H^^^ be- 
comes symmetric in this domain. The adjoint Hamil- 
tonian has the same differential form ([3]), but the 
domain is given by (/'(I2I) S V{H*^^), where and 
are absolutely continuous. The two domains are 
seen to be different, i.e., ^ |). So we need 

to find out the deficiency space solutions (square inte- 
grable), which are in our case the solutions of the differ- 
ential equations ,(j)± = ±i(j)±. These are 



b±{\z\) = W^^^^{ ia+\z\)), 

Ka± 



(5) 



= {0(1^1) + + exp(zS)0_(|z|)} , (9) 

where (/id^l) G V{H\z\)- Now we can solve the eigenvalue 
problem ^ and find out eigenvalue using the extended 
domain ^ . The eigenfunction of ([3]) apart from normal- 
ization is given by 



x{\z\)^W^,^{ — {a+\z\))., 



(10) 



where a = 57Jy~"e~ ^^^d m has already been defined. 

The eigenvalue of (fTU)) can be found by looking at the 
short distance behavior of and PH)) and equating the 
coefficients of equal powers of (a-|- |z|). For general value 
of the self-adjoint extension parameter S, the eigenvalue 
equation is found to be 



_ r(l/2 + m - a) _ xi cos{0i - S/2) 



r(l/2 



a 



X2 cos( 



S/2) 



(11) 



where r(l/2- 
a+) = X2exp( 



TO — a^) = xi exp(— i^i) and r(l/2 + TO — 
—i02)- The eigenvalue equation ITTI) is a 
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1-parameter family of equation. It is to be noted that 
the equation of the form pT|) is discussed in detail in 
Ref. [l^, [2^ and the reason we get similar equation here 
is that the nature of the boundary condition obtained 
from the short distance behavior of the special functions 
are similar in both cases. Each value of the parameter S 
correspond to a different quantization for the system. A 
plot of J{E\z\) as a function of a can be found in FIG. 1. 
The intersection of the horizontal line with the curves be- 
tween consecutive vertical lines gives the value of alpha, 
in other words value of the energy -B|z|- It is in general 
hard to solve pT|) analytically. But in certain limits it 
can be evaluated exactly. For example, when the R.H.S 
of (fTTj) is zero, the eigenvalue is 



1 



4^2(1/2 - rn + nY 



(12) 



where n — 0,1,2,3,.... When the R.H.S is infinity, the 
eigenvalue E^^^ becomes 



where a = 

2k 



and C{ia) and D{ia) are undeter- 
mined coefficients. One can now determine the coeffi- 
cients C and D by looking at the short distance behavior 
of ^ and and equating the same powers of (a-f |z|), 
which will now depend upon the self-adjoint extension 
parameter E. Note that the scattering state solution (fTS]) 
does not belong to the domain ^ for the simple reason 
that the scattering state solutions are not normalizable 
and thus does not belong to the Hilbert space. But there 
is no problem for the scattering states to be square in- 
tegrable at short distance and thus it may belong to the 
domain © as far as the short distance behavior is con- 
cerned (see [ll| for detail discussion with an example) It 
is also possible to calculate the S'-matrix and phase shift 
from this scattering state solution. 



E, 



4^2(1/2 + m + nY 



(13) 



In between zero and infinity the eigenvalue can be ob- 
tained from FIG. 1. Let us now consider the infinite 
magnetic field limit, B ^ 00. In this limit m ^ 1/2 and 
from we get 



, n-0,1,2,. 



(14) 



Note that for B ~> 00, we get the same energy eigenvalue 
(fT4)l from Eq. (US]), but this time n = 1, 2, 3, ... 



V. SCATTERING STATE SOLUTION OF 
EXCITON 

We now move to the discussion of scattering state solu- 
tions. The eigenvalue equation for scattering states will 
be the same Eq. where E^^i is now positive. The 
scattering state solution is of the form 

Xd^l) = C{ia)Ahs.,m ( -^(a + k|) ) + 



D{ia)W,s.,m[-^{a+\z\)] (15) 
iKa 



VI. DISCUSSION 



The knowledge of quantum mechanical behavior of ex- 
citon in a magnetic field background is essential for the 
study of absorption in semiconductors as has been shown 
in Ref. [13] ■ Since the time independent Schrodinger 
equation for the exciton cannot be exactly solved, an ap- 
proximation scheme is generally used in order to solve 
the system analytically. In order to overcome the short- 
comings of the approximation scheme and to guarantee 
unitary evolution of the system we obtain a 1-parameter 
family of self-adjoint extensions of the effective one di- 
mensional eigenvalue equation. We obtain bound state 
and scattering state solutions, which in general depend 
on the self-adjoint extension parameter E. The direct 
optical absorption coefficients K{T,) depend upon the 
bound and scattering state solutions which are depen- 
dent on the boundary conditions characterized by the 
parameter E. Therefore the parameter S can be used 
to fine tune the optical absorption coefficient K{Yj) to 
match with the experimental data. Our procedure gives 
all possible solutions for the system. 
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